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Abstract 

The AdS/QCD models are known to be tightly related with the 
QCD sum rules in the large--/V c (called also planar) limit. Rewriting 
the theory of infinite tower of free stable mesons expected in the large- 
N c QCD as a five-dimensional theory we scrutinize to what extend the 
bottom-up holographic models may be viewed as an alternative lan- 
guage expressing the phenomenology of planar QCD sum rules. It 
is found that many features of AdS/QCD models can be thereby ob- 
tained without invoking prescriptions from the original AdS/CFT cor- 
respondence. Under some assumptions, all possibilities leading to sim- 
ple Regge trajectories are classified and it is argued that the most phe- 
nomenologically consistent model is the one called "soft wall model" 
in the holographic approach, with a preference to the positive-sign 
dilaton background. 



1 Motivation and formulation of problem 



In the last several years, numerous attempts to apply the ideas of AdS/CFT 
correspondence [HE] to QCD have formed a rapidly developing line of re- 
searches in the modern field theory and phenomenology The main efforts 
are concentrated on description of the hadron and glueball spectra, the chi- 
ral symmetry breaking, and the related phenomenology [5H4"2"]. although the 
range of applications is now much broader. The holographic methods are 
interesting as a language that allows to discuss within a uniform framework 
various approaches to modeling the interactions and spectral characteristics 
of light hadrons, heavy-light systems, hadron form factors, QCD phase dia- 
gram, and other phenomenological aspects which were previously the subjects 
of investigations for different communities. This unifying property brings 
up a natural question — why the AdS/QCD models are rather successful 
at reproducing the non-perturbative physics of strong interactions and how 
the holographic language is interrelated with other known phenomenological 
methods? At present there are plenty of papers devoted to the applications 
of holographic ideas but there is a lack of works which would try to answer 
this question!]]. The given research is intended to shed light on the problem 
by observing that rewriting the large-iV c QCD in a fixed channel as a free 
five- dimensional theory, one is able to obtain models which are very similar 
to the holographic models of the bottom-up approach. 

The AdS/QCD models are essentially based on the identification of 4D 
Kaluza-Klein (KK) harmonics of a 5D field with an infinite tower of colorless 
free meson states that is expected in the large- N c (called also planar) limit of 
QCD [43HH]. Such an identification can represent a convenient mathematical 
trick: Any field <p(x, z) living in a compact 5D space can be decomposed in 
4D harmonics, (p(x,z) = J2^=o V 9 «(^)V'n(^), which are 4D fields with masses 
determined by topology of the 5D space. For this reason we may rewrite 
an infinite tower of 4D fields ip n {x) with equal quantum numbers and with 
a given spectrum m n as a single 5D field <p(x, z) propagating in a suitable 
background. Decomposing this field in 4D harmonics with normalizable co- 
efficient functions f n { z ), the dependence on the fifth coordinate z can be 
integrated out and one comes back to the original action for the tower of 4D 
fields ip n {x). The AdS/CFT conjecture, however, has much reacher dynam- 
ical content than just KK- rewriting — it asserts an exact correspondence 
between the generating functional of correlators in the 4D theory and the 

1 A likely exception is the paper [42] where it was demonstrated numerically that the 
phenomenology of spin-1 mesons within the hard wall holographic models depends slightly 
on the choice of metric, the AdS metric, however, minimizes deviations from the experi- 
mental data. 
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effective action of the 5D theory in which the ultraviolet boundary values 
of the 5D fields are identified with the sources in 4D theory [2J. The ques- 
tion is whether we may interpret the KK-rewriting in such a way that it 
mimicked the correspondence above? The answer seems to be positive as 
is demonstrated in Section 2. The idea is that after integration over z one 
obtains also boundary terms which may be identified with sources coupled to 
the fields ip n (x). This identification yields the same results for the two-point 
correlators as within the holographic models. 

In Section 3, under some assumptions we analyze all possibilities for the 
5D backgrounds which lead to simple Regge like spectra and then, in Sec- 
tion 4, we select the most reasonable 5D rewriting of large- iV c QCD basing 
on requirements to reproduce the analytic structure of two-point correlators 
and the expected phenomenological form of Regge spectrum, this rewriting 
turns out to be nothing but the so-called soft wall model [H]. 

The description of the chiral symmetry breaking is among the most dis- 
cussed issues in the holographic approach [8j[9l[20l[2Tl[3l]. We propose in 
Section 5 a modification of method of Refs. [SJ[H] that seems to allow for 
embedding this phenomenon into our approach. 

In Section 6, we consider interrelation between holographic bottom-up 
approach and the phenomenology of QCD sum rules in the large- N c limit. 

In the concluding Section 7, we discuss briefly the relation of the 5D 
formalism to some traditional phenomenological approaches. 

2 Derivation of holographic like models via 
Kaluza-Klein reduction 

With respect to the number of colors N c , the meson masses behave as m ~ N® 
and their full decay widths do as T ~ iV" 1 . Since the meson masses vary 
slightly with N c , it is thus often useful in the problems of finding the meson 
spectrum to take the limit of large N c . In the extreme case N c — > oo, the 
mesons are infinitely narrow and non-interacting, in addition, their appear 
in infinite tower [131144] at fixed quantum numbers. Such an infinite tower 
of resonance poles saturates completely the two-point correlation function 
of quark currents with quantum numbers corresponding to the given tower 
(see, e.g., the representation (11061) in the vector channel). With the help of 
functional integral formalism, this situation can be formally described by the 
following action, 

/oo 
n=0 
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where we have included the couplings to external sources which appear in 
the corresponding generating functional of the connected correlators. Here 

0j = <f>nin3...nj, Hi — 0,1,2,3 with the signature (H ), corresponds to 

a meson field of spin J whose quantum numbers J, G, P, C are not specified, 
and Of is a source which can be represented as 

Of = F (J) 0jy (2 ) 
where are the decay constants defined by 

(0|OS B) |^ B) > = ^ (3) 

for a meson (jyf 1 with "polarization" Sj and Oj is a common source to which 
the states 0^ are coupled with a "coupling" Fn \ In case of high spin fields, 
J > 1, under <9 M we mean the (general coordinate) covariant derivative V M . 
Another subtlety is that for J > 1 mesons many additional quadratic terms 
appear in the action (III) , however, the auxiliary conditions can be chosen such 
that they do not contribute (to be discussed below). The representation (J2j) 
follows from the requirement to obtain finally the standard form for two-point 
correlators, namely integrating formally over the field <j)j in the generating 
functional 

Z[Oj]= [ D<Pje T ^, (4) 



and differentiating twice with respect to Oj at Oj = one arrives at the 
sum over meson poles, 

(Oj{q)Oj(-q)) ~ 2 \ — V ( 5 ) 
^ q z — m„ T + is 

n=0 n ' J 

for the two-point correlators in the momentum representation. As is usual 
in description of higher spin massive bosons, the tensor <pj is symmetric, 
traceless, 4>^ = 0, to provide the irreducible (| J, \ J) representation of the 
homogeneous Lorentz group (or, equivalently, the positivity of energy), and 
satisfy the auxiliary condition <9 M M ... = to give the required 2J + 1 physical 
degrees of freedom. The tensor structure in the r.h.s. of representation (jSJ) 
depends on the structure of terms with derivatives in the action (JTJ), it will 
be not of interest for us. 

We will not consider the baryons because their masses behave as m ~ N c 
and their are not narrow at large- N c . The baryons emerge likely as solitonic 
objects in the large- N c limit [H], our further discussions are definitely not 
applicable to such a case. 
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Our task is to rewrite the action ([T]) as an action of some free 5D theory 
that is not necessarily covariant with respect to the fifth coordinate z, i.e. in 
the form 

S 5 d = (-1) J / d A xdzh(z) {d MVj d M y J -m 2 jf 2 (z) V j V J ) , (6) 

where M = 0,1,2,3,4, cpj = <fj(x,z), and fi(z), f2{z) are yet unknown 
functions of the fifth (space-like) coordinate z. We require that the action (jH]) 
must respect the invariance under the 5D general coordinate transformations 
at vanishing \z\. This implies two principal consequences. First, 

h{z) = e*My/\detG MN \ = e*Ma 5 (z), (7) 

where Gmn is the metric tensor defined from the 5D metric, 

ds 2 = a 2 (z)(dx ll dx fJl - dz 2 ) } (8) 

and possible deviations from the 5D covariance at large enough \z\ are pa- 
rametrized by means of a general factor e*'- 2 -', such a parametrization is 
convenient to draw parallels with the soft wall AdS/QCD models. Second, 
one has to respect the rule for contraction of indices, for instance, 

d M yjd M y J = d M VM 1 ...M J dM>VM>...M' J G MMl G M ^ . . . G MjM ' J , (9) 

where 

G MN = G- M \ = a- 2 (z) v MN . (10) 

We remind the reader again that Dm denotes the covariant derivative Vm for 
high spin fields. Finally, the action is (up to a general normalization factor 
which is set to 1) 

S 5 d = (-1) J | d 4 xdze^a- 2J+3 (z) {{d^jf - (d z vj) 2 -m 2 a 2 (z)ri} . 

(11) 

We will regard the 5D fields with J > as the gauge ones and use this gauge 
freedom to go to such a gauge that no additional quadratic terms appear in 
description of free higher-spin mesons. For instance, in the case of AdS space 
(or at least the AdS space in the UV limit), this is the axial gauge [i~H l4"5]. 

<Pz... = 0, (12) 

which will be implied in what follows. Together with the traceless condition 
this gauge fixes components of fields <fj. 
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Consider now the following Sturm-Liouville (SL) problem 

- d z [gj{z)d z ^\z)\ + gj {z)a\z)m^\z) = (13) 

with 

gj(z) = e^a- 2J+ \z). (14) 

We remind the reader the general results of SL theory in Appendix A. The 
equation (TT31) is the equation of motion for the action fTTTj) in the gauge ffT2]) 
if we admit 

<p(x, z) = e iqnX <p n (z), ql = ml. (15) 

The relation (fT5|) is usually interpreted as the string excitation corresponding 
to 4D particle with physical momentum q n . We will not use such an inter- 
pretation and instead of this assumption consider a general mathematical 
problem: Given a spectrum m? n j in the action ([T]), find the functions gj{z) 
and a(z) which being inserted in (1131) would give m? n j as eigenvalues of some 
SL problem. This constitutes the main intermediate step in rewriting the 
action ([1]) as some free 5D field theory. 

Under the conditions formulated in Appendix A, the SL problem (fl~3j) has 
the solutions Lfn\z) which are normalized as follows 



-max 



gj(z)<p£>(z)<pip(z)dz = 6 mn , (16) 

and form a complete set of functions, hence, the function ipj(x,z) in the 
action (fTTT) can be expanded in the 4D harmonics, 

oo 

Vj{x,z) = Y,^(x)^\z). (17) 

n=0 

Let us substitute the expansion (TT7I) in the action (TTT]) . 

/oo 
d±xdzgj(z) J2 {^^d^fd^f 

m,n=0 

where the notation f[T4l) has been used. Integrating by parts and making use 
of Eq. f[T3"l) . the second term in the action ([TBI) can be rewritten as (omitting 
the general factor — 0j m Vj^) 

max 

dzgj{z)d M <p$d M <pW = 

iin 

/"■^m ax 

<p$gj{z)d,<p<p + / dvp$<pWgj{z) K, - a 2 (z)m 2 j) . (19) 
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Now we can integrate over z in the action (IT51) with the help of (flBj) and ffT^]) , 
the result is 



/DC 
n=0 

-^jW^jC^z)!*^}, (20) 



with ipj(x,z) given by (TTTf) . The action (1201) has the form of the action ([I]) 
if we identify 

Of = -gj{z)d z ^ Vj {x : z)\l—. (21) 

If = then we have just rewritten the 4D action through the 5D one as 
in the Kaluza-Klein reduction. But we want to have an exact correspondence 
between the 4D action (BQ) and the 5D action (TlTi) . i.e. O^f 1 ^ 0. By choosing 
/3 = 0or/3 = 7r/2in the boundary condition (A. 3) (see Appendix A) we 
can nullify the term at z = z max in the Eq. (l2~Tj) . Note in passing that for 
the vector mesons the first possibility is actually realized in the soft wall 
holographic models [13] , where z max — > oo, while the second one — in the 
hard wall models [814T0]. where z max is the infrared cutoff. Thus, 



{ ?= lim gj (z)d,<piP(z)<pj{x,z). (22) 

Now we should identify the general source Oj such that the decay con- 
stants in the representation were non-zero finite numbers. In general, 
this identification depends on the choice of metric (the function a(z)) and of 
five-dimensional mass m 2 . If we fix the AdS metric at z — > z m m, o,{z) ~ 1/z, 
z > 0, and m 2 j in the form (1721) then we should identify 

Oj = zgjifj/c ~ JV-Vtpj/c, (23) 

F (J)= i im C ^^M ; (24) 

since the function (fif\z) in the realistic models (Eqs. (jUJ) and f loT]) ) behave 
as z 2 at small 2. The factor c appears because there is always the freedom to 
redefine F n — > cF n , O — >■ (9/c where the parameter c can be fixed by matching 
to the high energy asymptotics of the corresponding QCD correlators [46J. 
The expression (T241) reproduces the result of Ref. [8] for the vector mesons 
(up to a general factor depending on normalization of fields). 

Let us now compare the Kaluza-Klein like calculations above with the 
technique used in the standard holographic approach. Within the latter, one 
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evaluates the action ffTTj) on the solution of Eq. (TT3T) that would give in our 
notations 

S 5 b = (-1) J J d 4 xgj(z)d zV J Vj \ z=Zmia . (25) 

Supposing that ipQ is the Fourier transform of the source of QCD operator O j 
at the UV boundary and accepting <p J (q,z) = ip(q, z)(fQ (q) with the bound- 
ary condition 2 J ~V(<?> z)|*-h-o = 1 (we generalize the matter to the case of 
arbitrary spin J) one differentiates the Eq. (1251) twice with respect to the 
source and arrives at the same sum over meson poles as in the expression (JH]) 
using the spectral expansion of the Green's function for the Eq. (fT3]) , 

^ q z — mf lJ + %e 



n=0 



This very coincidence enables the assumed holographic duality j2] 

Z 4D [Oj] = S 5D [<p J (x,z min )] at z^ip J (x, z min ) = Oj (27) 

at the level of two-point correlation functions in the approximation of classical 
field theory. 

It should be noted finally that we did not need the interpretation of 
the fifth coordinate z as inverse energy scale that is used in the standard 
holographic approach. 



3 Regge spectrum 

On the theoretical and phenomenological [37] grounds, the spectrum of light 
mesons is expected to be of the Regge form, 

m 2 n j ~An + BJ + C, (28) 

where A, B, C are parameters, J is the spin, and n is the "radial" quantum 
number that enumerates the states with identical quantum numbers lying on 
the higher daughter trajectories. In addition, one expects A B, i.e. the 
behavior (see, e.g., (H1H9] f° r the recent discussions) 

<j ~ n + J. (29) 

In the given Section, we examine systematically at which conditions the 
Regge spectrum can be obtained. 
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For further analysis it is convenient to cast the Eq. (fT3"|) into the form of 
the Schrddinger equation via the substitution (we omit below the index (J) 
and the notation of dependence on z) 



ip n = e"*/V- 3 /Vn. (30) 



The result is 



-il% + Uif) n = m n if) n , (31) 
u" + (*)' + (* j) W + "" + a- J )^ +aX , ,32) 



2 V 2 J \2 

where the prime denotes the derivative with respect to z. As is known from 
Quantum Mechanics, the linear dependence rn? n ~ n takes place if U ~ z 2 at 
least at large z, i.e. if one has the oscillator type of potential. On the other 
hand, the Regge behavior 3 ~ J holds if there is a shift of energy levels 
(= masses square) rr? n ~ n + C(J) such that C(J) ~ J. With the help of 
formulae in Appendix B, the both conditions can be achieved simultaneously 
at many choices of the functions &{z) and a(z) in the potential (132]) . This 
freedom can be restricted significantly if we make the following simplifying 
assumptions: 

(a) The functions $ and a, i.e. the shape of dilaton and the metric in the 
AdS/CFT language, do not depend on spirH J; 

(b) The functions and a are continuous and differentiate on the interval 
(0, oo) or (— oo, oo) (depending on a model); 

(c) Subleading in J and n corrections to the Regge spectrum are absentH. 

Under these assumptions we are able to classify all possible models yield- 
ing the Regge like spectrum. 

The most important class of models that we will refer to as models I 
corresponds to the choice, 

I: $ = ±X 2 z 2 ; a = (R/z) k , where < z < oo. 
The parameters A and R have dimension of mass and inverse mass corre- 
spondingly. Without loss of generality, we set A = 1 and R — 1. The 
dependence on A and R can be restored in the final expressions by replacing 



2 Although there are speculations in the literature on that point |27j . 

3 In reality, this seems to be not the case. The form of these corrections, however, is 
not known even approximately, it depends strongly on a model and is highly speculative 
in the phenomenology. We would not like to boil down into such speculations. 
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z — > Xz, rrij — > mjR k , m n — > m n /X. The equation f[3~Tl) - f[3"2"|) reads then as 
follows 



2 2 



-<+^ 2 ±[^(2j-3) + i]+H^-JLi] — * + 



mj) n . (33) 



The normalizable solutions of Eq. (133]) under different subcases are given 
below (see Appendix B). 

IA+: $ = z 2 ; a = z~ k , k > 1; mj = 0. 
The spectrum is 

m 2 = 4(n + 1) + 2(|e M | + £ k ,j), n = 0, 1, 2, . . . , (34) 

where 

£m = *(-J-![)-^ ( 35 ) 
The eigenfunctions (see the notation f[3"U]) ) 



^^/^b e "^' hilJ ^ Jl ^- (36) 



The spectrum is absent for 



1< J <14 < 37 > 



The restriction ( 157|) means that at 1 < k < 2 there is no finite discrete 
spectrum for J = 2 states. The case k = 2 is special, the equation f l3"3"|) is 
then 

-<+(z 2 + 3)^ = m^n, (38) 

that yields formally the spectrum of masses and eigenfunctions given by 
Eqs. (B.4) and (B.5). We must remember, however, that, first, due to the 
pole in the metric the problem is defined in < z < oo while the eigen- 
functions (B.5) are normalized in — oo < z < oo, second, we should have 
<^n(0) = 0, i.e. only odd n in (B.5) should be selected. Keepeng this in mind, 
we arrive at the spectrum 

J,k = 2: m 2 m =4(/ + 3), 1 = 0,1,2,..., (39) 

and normalized eigenfunctions 



J, k = 2 : cp 2l+1 = ^J- ( -^ e ^z- 1 H 2l+l {z). (40) 
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IA : $ = -z 2 ; a = z~ k , k > 1; rrij — 0. 
The equation f l3"B"j) yields the spectrum 

m 2 n = 4(n + 1) + 2(|£ M | - n = 0, 1, 2, . . . , (41) 

and the eigenfunctions 



^'(dfW^ 1 ^- (42) 

As before, the spectrum of J = 2 mesons is absent for 1 < k < 2 and for 
k = 2 is given by 

J,fc = 2: m 2 m =4Z, Z = 0,l,2,..., (43) 



7T- 1 / 4 / 1 

J, fc = 2 : = —J —-^ z ^i(z). (44) 

IB+: $ = z 2 ; a = l/z. 
The equation f l3"3"|) is 

- < + + 2(J - 1) + (J " 2)8 + m '" 1/4 }^„ = m^n- (45) 
The spectrum 

m 2 =4n + 2(0 + J), ra = 0, 1,2, . . . , (46) 

where 



0= X /(J-2) 2 + m 2 J . (47) 

The eigenfunctions 



The spectrum is absent for < £j < 1/2 and represents a special case for 
£j = 1/2 that can be analyzed as in the model IA for any concrete value of 

2 

rrij. 

IB : $ = -z 2 ; a = l/z. 
The equation (1531 is 



< + {^- 2(J - 1)+ ^-^^-V4 | 



V>« = ™ 2 Vv (49) 
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The spectrum 

m 2 n = 4(n + 1) + 2(0 -J), n = 0, 1, 2, . . . , (50) 

= \Z3S z£j+Wl » (A (51) 

As before, the spectrum is absent for < £j < 1/2 and is special for £j = 1/2. 

IC ± : $ = ±z 2 ; a = z~ k , < k < 1. 
For mj = 0, this model is identical to the IA^ one but with one important 
distinction: The interval of forbidden spins (|37|) is now larger, in the limit 
k +0 only the scalar and vector modes can be described. In the case 
mj 7^ 0, we are not aware of analytical solutions, presumably they violate 
our assumption (c). It is clear, however, that at small enough and large 
enough z the solutions are approximately given by the mj = case. 

ID ± : $ = ±z 2 ; a = z~ k , jfe > 1; mj ^ 0. 
The spectrum is given by 

m 2 n = 4n±fc(2J-3) + l + ..., n = 0,1,2,..., (52) 

where the corrections presumably violate our assumption (c). The eigenfunc- 
tions are not known, at large enough z they are approximately given by the 
model IA ± . 

The second class of models that we will denote the models II corresponds 
formally to the limit k = in the models I, namely 

II: $ = ±X 2 z 2 ; a = const, where — oo < z < oo. 
Without loss of generality, we may set A = 1 and a — 1. Note that the 
absolute value of m 2 in the potential (I3"2"|) is not fixed because of freedom in 
choosing a. The model II has two different variants which are given below. 

II+: $ = z 2 . 
The spectrum is 

m 2 n = 2(n + 1) +m], n = 0,1,2,..., (53) 
i>n = ^f^-V'e-'H^z). (54) 

II": $ = -z 2 . 
The spectrum is 

m 2 n = 2n + m 2 j, n = 0,1,2,..., (55) 
^ = y^~ 1/4 ^W- (56) 



12 



There is the third, rather exotic class of models, it would correspond to 
the choice k = — 1 in the model I but with a broader possibility for the choice 
of function $, 

III: a = z/R. 

The models III can be divided into the following subcases. 

IIIA ± : $ = ±AV; A 4 + m 2 j/R 2 > 0. 
The equation (I3TI) is 

- € + { (A 4 + |f ) z 2 ± 2A 2 (2 - J) + {J ~ l? z ~ l,A } Vn = mjty„. (57) 

The slope of Regge trajectories includes the 5D mass mj, hence, this quantity 
should not depend on J. Without loss of generality, we may set X^+m^/R 2 = 
1 and R = 1. The equation (1571) does not describe the vector mesons, for 
J ^ 1 the spectrum is 

m 2 = 4n + 2| J - l| + 2 ± 2A 2 (2 - J), n = 0,1,2,..., (58) 



An interesting particular case here is the choice A = that corresponds to 
the absence of dilaton background in the holographic models, the slope of 
Regge trajectories is then determined by the universal 5D mass square m 2 . 

IIIB: $ = 61og(A^); (mj/R) 2 > 0. 
Choosing the units A = 1, the equation fl3Tl) reads 

1 I a ~\ 

t/j n = m 2 J) n . (60) 



- < + {(^)V + ^-i-*>;+>-V4 } 



Here the slope is always determined by m 2 . Setting for simplicity (mj/R) 2 = 
1, the spectrum is 

m 2 n = An + 20,b + 2, n = 0,1,2,..., (61) 

where 

6,6 = v/(J-l-6) 2 + 6, (62) 
and the eigenfunctions are given by 



e- z2 / 2 ^- 1+ ^^ 6 / 2 L^(z 2 ). (63) 



This model describes the mesons of all spins if b > (y/2- l)/2; at b = 
(\/2 — l)/2, the spectrum of J = 1 mesons is given by Eqs. (B.4) and (B.5) 
with the same reservation as for the model IA. 
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4 Choosing the most viable model 



After having classified all 5D models leading to simple Regge trajectories 
we should impose some criteria which would select the most viable model. 
Our first criterion is that a model should describe the mesons of all spins in 
a uniform way, i.e. without "pathologies" at certain values of J. On this 
ground, we reject the models of type IIIA and of type IC. 

The models in question can be matched to QCD by means of the quark- 
hadron duality principle that is considered in Section 6. The requirement of 
such a duality for the vector mesons constitutes our second selection criterion. 
This kind of matching for the present models cannot be fulfilled for the 
J > 1 mesons (at least we could not find the relevant examples), this is not, 
however, a serious drawback since the issue of quark-hadron duality for the 
J / 1 mesons is questionable by itselfl Let us show that the models of 
type II do not pass this criterium. 

Within the models under consideration, the condition of duality (IllOp 
translates into 

lim [d z ip n (z)f ~ m 2 n d n m 2 n} n > 1, (64) 

z— s-0 

where the l.h.s. is the residue in the case of type II models and (p n (z) ~ H n (z) 
(see Eqs. (|54|) and (|56p ). First of all, in order to have non-zero residues, 
we must define the type II models on the interval < z < oo, this changes 
slightly the normalization of wave functions (|54|) and (|56|) . Due to the identity 

d z H n (z) =2n# n _i(2), (65) 

only the states with uneven n survive: n = 21 + 1, I = 0, 1, 2, . . . . Making 
use of the property 

H 2l {0) = {-l) l 2 l (2l-l)\\, (66) 
and the Regge form of the spectrum, ~ n, we arrive at the expression 

[^ 2m (0)f ~ ml l+1 h h h = m ~^ ll]2 ■ (67) 

The quantity hi is a degreasing function of index / while it must be a constant 
(at least for large I) for the Regge like spectrum. Thus, the condition of 
duality (|61|) is not fulfilled. 

The condition f)64p is also very restrictive for the models of type I. First 
of all, it is easy to see that the residues are finite and non-zero constants only 

4 In the sense that the partonic logarithm given by the free quark loop does not yield a 
clear-cut main contribution to the two-point functions. 
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if the 5D vector field is massless. In this case (see Eq. (j55|) and expressions 
for the relevant wave functions) 




ymg(z)d,<p n (z) ~ \lj^m Ln2 (0) ~ V J^WY V » 



-. (68) 



nl 

Substituting this expression in the duality condition fl64|) we conclude that the 
latter holds for k — 1 only, i.e. for the AdS background in the limit z — » +0. 
This situation corresponds to the models of type IB with rrij =1 = 0. 

Our third selection criterion is the requirement to have in a natural way 
the spectrum ( 1291) . i.e. the slopes of spin and radial trajectories must coincide. 
We should emphasize a remarkable fact that the third and the second criteria 
have a significant overlap. For instance, the form of spectrum ( 129]) selects 
immediately the models IB among the models of type I. Also it excludes the 
model IIIB. Thus, the phenomenological spectrum (129 p and the quark-hadron 
duality principle strongly suggest that the AdS background in the limit z — > +0 
is the only possible 5D realization of large-N c QCD. 

We note further that the spectrum depends on spin J not only by con- 
struction (due to the covariant way for contraction of Lorentz indices) but 
also through the dependence on J of 5D mass rrtj. The latter dependence 
is ad hoc and requires invoking some speculative assumptions. It would be 
natural therefore to have a minimal impact of this dependence on the final 
results. In the models of type II, the spectrum depends on J completely 
through m 2 , such a realization of Regge trajectories looks unfortunate and 
for this reason it does not worth serious consideration. 

The third criterion might be suggestive in the choice between the mod- 
els IB + and IB . In order to have a simple Regge-like spectrum (our as- 
sumption (c)) 

mlj = 4(n + J + C), (69) 

where C is a constant and the value of the slope was determined by the 
choice A = 1 in the dilaton background e ±A z , we must impose the following 
dependence of 5D mass on spin J (see Eqs. ( )46|) . (|47|) . and ( 150]) ) 

IB+: m 2 j = 4J(C+1) + 4(C 2 -1), (70) 
IB : m 2 j = 8J 2 + 4 J(3C - 2) + 4C(C -2). (71) 

The condition m 2 J=l = obtained above yields for the both cases one ac- 
ceptable solution C = that is, however, not valid for the J = case in the 



15 



model IB . The five dimensional mass is then 

IB+ : m 2 j = A(J-l), (72) 
IB": m 2 j = 8J(J-l), J > 0. (73) 

The both expression result in the spectrum 

m 2 n>J = 4(n + J), J>0. (74) 

The spectrum (I74p turns out to be not valid for the scalars of the model IB + 
as well because the equation (143]) does not have finite discrete spectrum in 
this case, i.e. when m 2 J=0 = —4 according to Eq. ( 172]) . Thus the scalar 
mesons represent a special case in the models IB and should be considered 
separately. It seems to be natural to supplement the relations (172]) and ( 173]) 
by the condition^. 

m 2 J=0 = 0, (75) 

The requirement ( 175]) entails a change of prescription ( 124]) for the residues, 
namely the latter must be replaced for the scalar case by 

F n = hm c — . (76) 

The condition (175]) leads to the following spectrum for the scalar mesons, 

IB+: r<o = 4(n + l), (77) 
IB": m£ >0 = 4(n + 2), (78) 

i.e. the scalars are degenerate with the vector mesons in the model IB + and 
with the J = 2 tensor mesons in the model IB - . 

Consider now the impact of 5D mass mj on the final results of models IB. 
For this purpose let us just set mj = 0, the spectra then are 

IB+: m* J<2 = 4(n + l), (79) 

<j>2 = 4(n + J -I); (80) 

IB - : r< J<2 = 4(n + 2- J), (81) 

<j> 2 = 4n. (82) 

We see that the impact of mj is rather slight in the model IB + and is deci- 
sive in the model IB - . According to our third criterion, the model IB + looks 



5 Due to accidental coincidence of this requirement with the Eq. (175]) in the J = case, 
the condition (|75[) may be imposed on the model IB + only. 
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therefore more attractive phenomenologically than the model IB - . This con- 
clusion, however, may turn out to be misleading if we attempt to match 
the present phenomenological approach to the theory of free higher spin 
fields propagating in AdSs background where the mass coefficient behaves as 
J 2 — J — 4 (see, e.g., [4"5|l50] and references therein), the related discussions 
are beyond the scope of the present paper. 

Let us compare our results with the assumptions used in the holographic 
models. According to the principle of AdS/CFT correspondence, each oper- 
ator Oj{x) in the 4D field theory corresponds to a field ipj(x,z) in the 5D 
bulk theory that has a mass mj determined via the relation [2] 

m 2 j = (Aj- J)(Aj + J-4), (83) 

where Aj is the canonical dimension of the operator Oj(x). One can con- 
struct two types of minimal interpolating operators which describe the mesons 
in QCD, 

Of 2 = q{ lh ) 1{lxl D^.--D lxj} q, (84) 

O t J = 3 = q(l,)D { , 1 D, 2 ---D, j} q, (85) 

where t denotes the twist, t = A — J. If we restrict ourselves by the operators 
with minimal twist (a custom first approximation in QCD) then Aj = J + 2, 
hence, according to Eq. (I83p . we arrive exactly at the relation (1721) ! The op- 
erators of the minimal twist, however, cannot interpolate the scalar mesons, 
in this case one uses the interpolator fl85|) that predicts A = 3, hence, 
mj =0 = — 3. This does not contradict to our analysis where the scalar case 
was also special. We would then obtain for the model IB + 

m 2 J=0 = -3: r< = 4(n + l/2). (86) 

The phenomenology should tell us which spectrum of scalars is preferable, f[TT|) 
or (I8"6"j) . If the spectrum (1771) is definitely better phenomenologically then, 
within the holographic approach, one should seemingly use the interpola- 
tor (1841) for the scalars also, i.e. make contraction of indices and interpolate 
the scalars by the twist-four operator Oj^ = qi^f^jl^D^q, thus A = 4 
and m 2 J=Q = according to Eq. (1HB1 . The matter is however complicated 
by the fact that, firstly, this interpolator is equivalent^ to Cj^q = rn q q{^i^)q 
due to the QCD equations of motion and thus it can be neglected in the 
chiral limit, m q = 0, secondly, one can construct a purely gluonic twist-four 
interpolator (9j= 4 = OLsG 2 ^ that makes difficult to distinguish between the 

6 Neglecting the axial anomaly in the pseudoscalar isoscalar channel that is suppressed 
in the large- N c limit. 
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quark-antiquark states and glueballs in the scalar isoscalar channel within 
the holographic approach. 

The assumption of the AdS background in the limit z — > that is ex- 
ploited in the AdS/QCD models matches perfectly our results as was re- 
marked above. 

The spectrum f)74p was first obtained in the holographic model of Ref. [13] . 
We note that the sign of dilaton background was opposite in that model 
(corresponding to the model IB~ in our classification), this would lead to a 
different spectrum in our approach which would coincide with the Eq. ( J74j) 
for the vector mesons only (compare the Eqs. fT46l) and fl50l) ). The reason 
of this disagreement lies in a different way of introduction of higher spin 
fields^ in [13] . According to the arguments of recent papers [571 158] based 
on realization of confinement in the string picture, the dilaton background 
should have positive sign in the soft wall AdS/QCD models, i.e. the most 
consistent models are those of type IB + in our classification. We have arrived 
at the same conclusion without using any string arguments. 

5 Chiral symmetry breaking 

The incorporation of the Chiral Symmetry Breaking (CSB) into 5D ap- 
proaches is often a questionable subject. Nevertheless we will try to outline 
a possible description of the CSB phenomenon that is compatible with our 
scheme of 5D holographic like rewriting. 

First of all, since we have no quarks, no microscopic model of CSB is 
possible, the best we may do is to describe the consequences of CSB on 
the hadron level. The first such consequence is the phenomenological fact 
that the masses of parity partners which seemingly belong to the same chiral 
multiplet are quite different. To reflect this effect one should introduce a 
mechanism for this mass splitting. Within the 5D framework, the chiral 
symmetry does not exist at all because there is no analogue for the matrix 
7 5 in five dimensions, the CSB can be therefore simulated only indirectly 
by means of somewhat different description of states with equal spin but 
opposite parity. We will restrict ourselves by the vector and scalar sectors 
only. 

The second consequence of CSB consists in the appearence of massless 
(in the chiral limit) pseudoscalar meson due to the Goldstone theorem. The 

7 We have introduced the higher spin fields in a direct way as in paper |36j . In the 
latter paper (and analogous ones), however, the hadron spin is subsequently devided into 
the orbital momentum of quark-antiquark pair and the proper quark spin while we do not 
use such a model assumption. 
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both consequences should be described with the help of a mechanism that 
explains, e.g., why the relation m n = is naturally related with m 2 ai > 2m 2 
(instead of m ai = m p as naively expected from the linear chiral symmetry). 
The simplest such a mechanism that is realized in the AdS/QCD models and 
that we will consider as well is borrowed from the low-energy effective field 
theories: One introduces a scalar field X that acquires a non-zero vacuum 
expectation value (v.e.v.) X (z) and is coupled to the axial- vector field 
through the covariant derivative [HE], 

Scsb = J d 4 x J dzy/)G^\e*W (\D M X\ 2 - m 2 x \X\ 2 ) , (87) 

where 

D M X = d M X - ig 5 A M X. (88) 

As usual, it is enough to restrict ourselves by the quadratic in fields part 
in the action (187)1 because the equations of motion can provide a non-zero 
v.e.v. already in this case if the bulk space is curved. This simplification 
turns out especially convenient for us since it is much easier to integrate 
over z in order to see the equivalent 4D effective theory. The equation of 
motion determining the v.e.v. Xq is nothing but the SL equation ffT3]) for 
the massless scalar particle. According to a recipe based on the AdS/CFT 
correspondence [51J, it must behave at z = as 

X (z)\ z ^o = C 1 z + C 2 z\ (89) 

where C\ is associated with the current quark mass, C\ ~ m q , and C% with 
the quark condensate, C 2 ~ (qq)- This interpretation implies a somehow well 
established correspondence of the model to QCD. As we do not have such 
a correspondence, it is more honest to say that the incorporation of the 4D 
massless scalar particle can be related to the spontaneous appearance of two 
order parameters with mass dimension one and three and this property may 
be exploited to mimic the CSB. 

It is easy to establish a general self-consistency condition for such a de- 
scription of the CSB. Substituting tp = X = z h into the Eq. (JTSJ) with 
a(z) = z~ k and zero r.h.s. (k > 0, h > 0) and retaining the leading terms at 
z —7- 0, one has 

m 2 xZ h-5k = h ( h _ 1 _ 3^-3^2, (90) 

that yields immediately k = 1, i.e. this design of CSB is possible at m 2 x ^ 
only in the AdS background at z — )■ 0. The second condition is m 2 x = h(h — 4) 
that both for h = 1 and for h = 3 gives m 2 x = —3 in agreement with the 
Eq. flHBl at A = 3. The case m 2 x = we will briefly analyze later. 



19 



The fact that the pion does not belong to the corresponding linear pseu- 
doscalar trajectory and should be considered separately is in agreement with 
the phenomenology. The solution Xq(z), however, turns out to be non- 
normalizable in the known bottom-up holographic models. This is a serious 
problem for our approach because we cannot integrate over z and see the 
underlying 4D image of the description. A way out may be the following. 
The equation of motion for X represents a second order linear differential 
equation that has two independent solutions, 

X Q (z) = C 1 X 1 (z) + C 2 X 2 (z). (91) 

The solution X\ (z) spoils the normalizability at z — > while (in the soft wall 
models) X 2 (z) does at z — > oo. We do not see any reasons why X (z) must 
have everywhere a continuous derivative. Taking this into account, we can 
construct the following normalizable solution, 

X (z) = C 2 X 2 (z)\ z < Z0 + CrX^z)^ (92) 

with the continuity condition 

C 1 X 1 (z ) = C 2 X 2 (z ). (93) 

Another restriction on the input parameters comes from the normalization (|16|) . 

Cl / g (z)X 2 2 (z)dz + C\ \ g (z)X 2 1 (z)dz = 1. (94) 
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The system of equations (19^1) . (I9~4"j) allows to exclude two of three parameters 
C\, C 2 , Zq. 

If zq is associated with the inverse energy scale then it is natural to 
interpret z^ 1 as the CSB scale, Acsb — 47r/ 7r ~ 1 -j- 1.2 GeV. The physics 
of strong interactions is known to be substantially different below Acsb and 
above Acsb, the effective emergence of this scale represents in fact the third 
consequence of the CSB in QCD. This consequence of the CSB seems to be 
not incorporated into the existing AdS/QCD models. 

Let us consider a concrete example — the model IB + that has been se- 
lected above as the most viable phenomenologically. The corresponding so- 
lution for Xq(z) is known [38J, 

X x {z) = ze~ z2 U(-l/2, 0; z 2 ), (95) 

X 2 (z) =z 3 e- z2 M(l/2,2;z 2 ), (96) 
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where U and M stay for the Kummer function U and the Kummer function 
M correspondingly. They have the asymptotics 

X^z) > X 1 {z) > z 2 e- z \ (97) 

Z— >0 2 — >00 

X 2 {z) >• z 3 , X 2 {z) >• 1/Vtt. (98) 

2— >-0 2— >00 

The continuity condition (|93|) takes the form 

dUi-1/2, 0; zl) = C 2 z 2 M(l/2, 2; z 2 ), (99) 
while the normalization yields 

C\ \ z 3 e~ z2 M 2 {l/2, 2; 2 2 )c/z + C? / z^e^C/^-l^, 0; 2 2 )c/z = 1. 

JO Jz 

(100) 

The solution of sum rules (I120p gives a remarkable relation for the slope 
of meson trajectories, a ~ A^gg. Identification of z with A^g B means then 
in our units for z that z = 1/2. The system of equations (19"9~1) . fllOOp yields 
numerically C\ ~ ±1.2, C 2 ~ ±3.5. 

The substitution of solution (1921) back into the action (1H7|) gives the vac- 
uum energy. As follows from the asymptotics (198]) . the vacuum energy density 
e vac will be divergent logarithmically at z = 0, 

£vac = 96*2 log 2! + const. (101) 

By subtracting the divergence in Eq. (11011) . we obtain the renormalized vac- 
uum energy density. In the case of numerical example above it is 

4; e c n > w -29.3| z < .5 + 3.2| 2>0 .5 = -26.1. (102) 

The negative value for e-^^ is in accord with expectations from the gluody- 
namics [52J. 

Comparing the prescription (1891 with the solution (1921) and the asymp- 
totics (|98|) . we would say in AdS/QCD language that the CSB is mimicked 
in the chiral limit (zero current quark mass). This is self-consistent as long 
as we have assumed initially the existence of exactly massless scalar particle. 

Let us consider the case of vector mesons. First of all we note that the 
trick above does not allow to introduce the massless vector particles as bound 
states — one can check that although non-normalizable massless solutions 
exist for the spin-one mesons, both of them are non-normalizable at z — > 0. 
Thus, we are save from the appearance of massless vectors^]. 



8 We would note that there is a related logical uneasiness in some AdS/QCD models: If 
one introduces a non-normalizable scalar field, why is this forbidden in the vector channel? 
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The axial- vector sector is a traditional place for modelling the CSB physics. 
We will not boil down to the model building of this kind (the most popular 
way is the introduction of longitudinal component for Am and of the pion 
field 7r through the parametrization X = Xoe l2n ) since in the chiral limit it 
does not seam to bring new interesting results. An important consequence 
of the term (1871) in the action is that the masses of axial- vector mesons get 
shifted due to coupling to the v.e.v. X through the covariant derivative (1551) . 
Within the model IB + , the new mass spectrum can be in principle calcu- 
lated via the substitution m 2 — >■ m 2 + gfX^z) into the Eq. (j45|) . Although 
the exact analytic solution is not known, it is clear that the masses will 
be enhanced because the new contribution to the "potential" of Shrodinger 
equation (|45|) is always positive. In addition, from the solution ( l§2"j) and the 
asymptotics (1971) . (1951) follows that the holographic potential in the Eq. (Ho^) 
will be (taking into account that m 2 J=1 = 0) 

U\ z » zo = z 2 (l + glCle- 2 * 2 ) + 1 (103) 

U\ z<<Z0 = z 2 (l + g 2 C 2 2 z 2 ) + ^. (104) 

This behavior demonstrates that the spectrum of axial-vector states ap- 
proaches rapidly to the spectrum of vector mesons and the axial-vector 
residues are well-defined at least for high enough excitations. The rate of 
such a "Chiral Symmetry Restoration" (CSR) is exponential that is in agree- 
ment with the analysis [531154] based on the QCD sum rules. This qualitative 
feature differs from the previous results obtained within the soft wall holo- 
graphic models [34"ll38] in which the shift between the V and A masses square 
tends to a constant and the axial- vector residues cannot be determined from 
the expression f[2~4"l) . It must be noted that the described behavior does not 
contradict to the CS non-restoration scenario observed in the phenomenol- 
ogy [SJ49J. The latter states that the mesons on the leading Regge tra- 
jectories do not have parity partners and the spectral degeneracy has the 
form 

m n+i,j~ m n,J U uneven); m+ +1J « m~j (J even), (105) 

where ± refer to parity. The model above can be accommodated to such a 
scenario. 
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6 QCD sum rules in the large-iV c limit vs. 
AdS/QCD 

The momentum representation for the vector (V) and axial-vector (A) two- 
point functions in the large-iV c limit of QCD can be written as (neglecting 
the nonpole terms) 

oo p 2 

<•#(?)•#(-?)> = KM) E C1 2 J'" F > ( 106 ) 

n=0 V m V,n <£ 



(Jtwft-q)) = + KM) E q2 m t w +ie . ( 107 ) 

9 n =o 5 A >™ 

where _F n are the corresponding electromagnetic meson decay constants, f n 
denotes the pion weak decay constant, and 

KM) = + Qj f- (108) 

is the transverse projector. 

On the other hand, one can write the Operator Product Expansion (OPE) 
for those two-point functions at large Euclidean momentum, Q 2 = —q 2 , 

( A 2 °° C \ 

(JM)U-q)) = KM)Q 2 [Co in— + E ^ J ■ ( 109 ) 

Here A is a renormalization scale and the coefficients C{ depend weakly on A 
and Q 2 . It is easy to show that in order to obtain the analytic behavior (I109p 
the residues in Eqs. (I106p and (I107p must satisfy the following condition (see, 
e.g., |53J for a detailed discussion) 



^2 



m 2 n d n m 2 n . (110) 



Thus, given the relation (IllOp . the infinite number of meson poles is equiv- 
alent to the partonic logarithm in the Eq. (11091) plus some corrections. This 
property we refer to as the quark-hadron duality in the large- N c limit (brief 
reviews can be found in [53TI55] ). Matching the Eq. f)109p to the represen- 
tations f)106p and (I107p in the Euclidean region, one obtains the QCD sum 
rules which represent certain relations on the parameters of the meson spec- 
trum P2H57J. 

Consider the simplest example — the linear spectrum 

mv A . n = a(n + b v ,A;n), (HI) 
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where the slope a is assumed to be universal. According to the condi- 
tion f lllOp . the residues must behave as (we neglect possible corrections to 
the relation ffTTU]) [53]) 



172 2 T72 

t V,A;n ~ m V,A\n*V,Ai 



F V,A 



const. (112) 
Performing the summation in n in the expression f)106p . one arrives at [51]loT] 

, (113) 



9 F 2 



a f (-1)V5^ 
7t2 + z__, h const 



i=l 



*g 



where Bi(x) are the Bernoulli polynomials. The OPE for the two-point 
functions at N c 3> 1 and in the chiral limit reads as follows [58] 



12tt 2 



7T 



A 2 



1 + — In— + 



(G 2 ) , 4£ 



Q 2 12tt g 4 



+ 



9 



-7ra, 



(gg) S 

g 6 



(114) 



Here (G 2 ) and (gg) denote the gluon and quark condensates and £ y = —7, 
£ A = 11. Comparing the Eq. (11141) with (II 13p and with analogous expression 
for the axial-vector two-point function, one obtains the following sum rules 
(we consider the leading order in a s ), 



F 2 =F 2 = F 2 







N c a 
24tt 2 


. 


(115) 


b v - 


1 

2' 


(116) 


b A - 


1 

2' 


(117) 


K,a- 


- K,A + -, 




(118) 


^v,a{ 


b V ,A - 1/2) (b VA - 1) . 


(119) 



f 2 

a s {G 2 } 
12vrF 2 a 

3F 2 a 2 

The sum rules fU15l) - flll9p cannot be satisfied by the simple linear spec- 
trum fillip , a natural solution of the problem is to consider nonlinear correc- 
tions to the ansatz fillip or/and not place some states on the linear trajectory. 
We note, however, that the sum rules ^J^ -^J^,^^- ^J3^ - ^J^ -. . . are 
progressively less and less reliable due to the asymptotic nature of OPE and 
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growing anomalous dimensions of condensate terms. The simple linear spec- 
trum is able to satisfy the most reliable relations (|115p - (|117j) . Identifying 
F = F p and making use of the KSFR relation [59] F 2 = 2f% as in the 
original Weinberg sum rules [60], we have immediately the solution 

b v = 1/2, b A = l, a = AM 2 fjN Cl (120) 

which is in a good agreement with the phenomenology (for instance, the 
Weinberg relation = 2m 2 p follows automatically). This solution for the 
intercept by a yields the spectrum of the generalized Lovelace-Shapiro (LS) 
dual amplitude |61J. The difference b A — by emerges from the CSB dynamics 
that is directly seen in the Eq. (1117P (namely the fact that f n ^ 0). Within 
the LS amplitude, the same difference arises from taking into account the 
Adler self-consistency condition [62] which is a logical consequence of the 
Goldstone theorem. 

Turning back to the holographic approach, we note that the spectrum ( 174)) 
that often appears in the soft-wall AdS/QCD models, in fact, does not de- 
scribe the p and u> mesons as opposed to the usual belief, it corresponds 
rather to the axial-vector mesons. We are not aware of any natural way for 
accommodation of the intercept by = 1/2 within the holographic approach 
and/or the technique presented in the previous sections, unless some specu- 
lative extra assumptions are involved [2T|[57] . As follows from the discussions 
above, the correct vector spectrum should naturally stem from a successful 
implementation of the CSB physics in the holographic model^E 

As seen from the expansion (11131) . the spectrum of masses (at least in 
the vector channels) determines the gauge-invariant QCD condensates; if 
there are several such condensates of equal dimension it determines a certain 
linear combination of them. For this reason, any attempts of inclusion of 
QCD condensates into the soft wall holographic models face a serious prob- 
lem of double counting. This remark, however, is not relevant for the hard 
wall models [8HTU]. where the masses are given by roots of Bessel function, 
m n ~ n, and for this reason it turns out that the corrections to the partonic 
logarithm does not reproduce the analytic structure of OPE. 

Our previous discussions incline us to conclude that the bottom-up holo- 
graphic approach may represent just an alternative language for expressing 
the phenomenology known from the QCD sum rules in the large- N c limit 
with practically the same number of input parameters and ensuing accuracy. 

9 The situation resembles the almost 40 years old discussions about a still unresolved 
problem: In contrast to the Veneziano dual amplitude that gave rise to the modern string 
theory, the LS amplitude does not ensue from any string model; the incorporation of the 
CSB effect into the string approach is a formidable task up to now. 
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This is directly seen from the bottom-up derivation of AdS/QCD like models 
in Section 2 and a posteriori from the type of questions which are addressed 
within the holographic models. First of all, the problem of finding the fittest 
spectrum in the sum rules is reformulated within the AdS / QCD models as the 
problem of finding an appropriate 5D background and boundary conditions 
on the 5D fields which give the required spectrum as a result of solution of 
the corresponding SL problem while in the sum rules the spectrum is looked 
for empirically ignoring its model-dependent origin. The subsequent phe- 
nomenology, however, is very similar in the both approaches. For instance, 
both sum rules and most of holographic models incorporate the pion via the 
Partial Conservation of Axial Current hypothesis. As is known from the phe- 
nomenology of sum rules [53], it is difficult to accommodate the pion as the 
first state in the tower of pseudoscalar resonances, one should rather consider 
it separately as a special state, the same conclusion seems to follow from the 
AdS/QCD models. 

7 Concluding remarks 

The phenomenology expected in the large- N c limit of QCD can be compactly 
rewritten in terms of some phenomenological 5D theory, namely the infinite 
towers of narrow meson states coupled to external sources can be formally 
represented as some 5D single fields propagating in a suitable background. 
This allows to replace sometimes lengthy manipulations with infinite series 
of meson poles by compact operations with those 5D fields. Needless to say 
that such an introduction of 5D formalism represents a certain operational 
progress and one could even draw a parallel with the introduction of con- 
tinual integral into the field theory in the sense that manipulations in the 
latter formalism are usually understood as a compact formal way to write 
lengthy operations with series of perturbation theory. The modern 5D model 
building, however, involves various speculative assumptions borrowed from 
the AdS/CFT correspondence which look often unfounded in their extrap- 
olation to real QCD. In the present work, we have explicitly demonstrated 
that 5D models very similar to the so-called soft wall holographic models 
can be derived within the 5D formalism imposing some requirements of phe- 
nomenological consistency without use of AdS / CFT dictionary. These mod- 
els represent just another language for expressing the Regge phenomenology 
within the framework of QCD sum rules. 

The 5D formalism proposes new ways for description of the chiral sym- 
metry breaking. The related attempts cannot be viewed as another language 
for expressing the matching of low-energy effective field theories for strong 
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interactions to QCD sum rules (see, e.g., [63] among others), but are similar 
in spirit though. The reason is that the effective models deal with the lowest 
mass (the first radial excitation at best) mesons while the 5D models are able 
to describe the impact of CSB on the whole tower of excited states. Another 
reason consists in the fact that the existing 5D descriptions of CSB are vague 
from the 4D point of view. We have tried to advance in "visualizing" a 5D 
description of CSB in 4D terms but could not provide a concrete example. 
The main problem to be solved is the derivation of GOR relation for the pion 
mass within the presented approach. 

We did not discuss an explicit incorporation of quarks and QCD running 
coupling since without knowledge of the confinement mechanism, at least in 
the planar limit of QCD, any such attempt is doomed to be rather specula- 
tive. 

It must be emphasized that the considered 5D approach is completely 
phenomenological and hardly can be directly related to any fundamental 
extra-dimensional theory. First of all, the 5D background of the soft wall 
model that we have also obtained cannot be a solution of 5D Einstein equa- 
tions as it breaks explicitly the 5D general covariance. A possible way out 
could consist in constructing an alternative mechanism for mass generation. 
Such an attempt for flat metric (an "effective" holographic model) was un- 
dertaken in Ref. [M] (see also [65]) where the Higgs mechanism was used. 
Second, we have ignored all long-standing problems in description of higher 
spin fields. If we included a gravitational dynamics that would yield particu- 
lar metrics in question we would have the problem of gravitational coupling 
for higher spin fields — such a theory does not exist. For this reason we 
believe that within the present approach it is meaningless to analyze the 
backreaction of metric due to condensation of scalar field when modeling the 
CSB physics as well as to address other questions of relation to consistent 
extra-dimensional theories until some fundamental problems in field theory 
are solved. Further discussions of relevant problems can be found in [66J. 
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Appendix A: Sturm-Liouville theorem 

In this Appendix, we briefly remind the reader the main results of the Sturm- 
Liouville (SL) theory. The SL equation is 

-d z [p(z)d z ip] + q(z)ip = \u(z)tp. (A.l) 

Here the function p{z) > has a continuous derivative, the functions q(z) > 
and ou(z) > are continuous on the finite closed interval [z m i n , z max \. Loosely 
speaking, the SL problem consists in finding the values of A for which there 
exists a non-trivial solution of Eq. (A.l) satisfying certain boundary condi- 
tions. Under the assumptions that p(z)' 1 , q(z), and u>(z) are real-valued 
integrable functions over the interval [z m [ n , z max ], with the boundary condi- 
tions of the form 

ip{z min ) cos a - p{z min )ip' '(z min ) sin a = 0, (A.2) 

y?(^max) cos (3 - p(z max )ip' (z mscK ) sin (3 = 0, (A.3) 

where a, (3 G [0, 7r) and prime denotes the derivative, the SL theorem states 
that 

• There is infinite discrete set of real eigenvalues A n , n — 0, 1, 2, 

• Up to a normalization constant, there is a unique eigenfunction <p n (z) 
corresponding to each eigenvalue A n and this eigenfunction has exactly 
n - 1 zeros in [z min , z mscx \. 

• The normalized eigenfunctions form an orthonormal basis 

ip m (z)ip n (z)u(z)dz = 5 mn . (A.4) 

- min 

Thus the solutions of the SL problem form a complete set of functions in the 
interval [z m i n , z max ], i.e. they can be used for expansion of arbitrary functions 
in that interval. 

Appendix B: Schrodinger equations appearing 
in the work 

The spectrum is typically determined by the equation 

i[) n = m 2 n ijj n , (B.l) 



< + 



z 2 + 



b 2 - 1/4 



+ c 
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that is exactly solvable. The spectrum of normalizable modes is 

m 2 n = 4n + 2\b\ + 2 + c, n = 0,1,2,..., (B.2) 
the corresponding normalized eigenfunctions are 



where l)n are associated Laguerre polynomials. 

In the case \b\ = 1/2, the solutions of Eq. (B.l) are different. The spec- 
trum is 

m 2 n = 2n+l + c, n = 0,1,2,..., (B.4) 
and the corresponding normalized eigenfunctions are given by 



if), 



where H n are Hermite polynomials. 

In the case |6| < 1/2, there is no finite discrete spectrum. 
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